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1. Introduction 

The 2D Ising model is certainly one of the most famous statistical models in physics 
with connections in various fields of research. Despite its apparent simplicity, it 
contains generic properties for phase transitions. From the interpretation of the Ising 
model as a lattice gas model, it was used as a simple model for wetting transitions. 
A number of theoretical and numerical tools have been developed to study such 
transitions. One of the first steps comes from McCoy and Wu who developed the 
theory of Toeplitz determinants Pj in order to solve the 2D Ising model with uniform 
boundary magnetic field (see also j^lElE]- Surface fields can be considered as the 
effect of boundary impurities or constraints that change the properties of bulk spins. 
They also can be viewed as a chemical potential difference between the bulk and the 
wall of the system in a binary mixture. The wall can indeed attract one species more 
than the other, leading so to an effective field 0. A major contribution is due to D.B. 
Abraham in the early 80 's who developed efficient methods based on transfer matrix 
techniques in order to solve the 2D Ising model with various fixed boundary conditions 
(the spins are constraint to be up or down on the boundary according to different 
profiles of the boundary spins). Theses techniques allow the computation of the surface 
tension energy for different profiles of the boundary conditions. Important questions 
arise whether the interface between two phases and produced by boundary impurities 
is diffuse or sharp (see discussion in 0) and how it diverges with the system size. These 
questions are closely related in general to the existence of roughening transitions in 
such systems close to the bulk critical temperature. A generalisation of McCoy and 
Wu result for the uniform boundary field case is to take two opposite surface fields 
Hi and H2 on a Ising strip, where a domain wall propagates in the middle of the strip 
when the fields have opposite signs [3|Hlini- A wetting temperature T^{Hi, H2) occurs 
at slightly below the bulk temperature where a well defined interface appears in the 
middle of the strip. This temperature can also be defined as the transition between 
diffuse and sharp interfacial regimes. In the particular case of an infinite strip, the 
width of the strip is only important for the magnetisation profile and for the scaling 
of the typical length of the interface extension [HlEB^DEl, which diverges typically 
like the inverse of (T — T^) above T^^, giving an interface correlation length exponent 
equal to unity. Technically, a generic way to implement a finite boundary magnetic 
field from an infinite one is to take the bonds perpendicular to the surface field line 
a fraction of the bulk bonds. For example if the spins subject to an infinite magnetic 
field Hi — ±00 and located along the surface arc noted (Ji^m n- = l.-.Ly, and if 
the bulk coupling is J, taking Jq as the coupling between (Ti,„ and (72,n generates on 
spins (72, Ti a field proportional to ±Jo/J. The value of this field can then be tuned by 
varying Jq , and this specificity is discussed in several publications |71 . In these 

different cases, the way the thermodynamic limit is taken is also important. 
The effect of an inhomogeneous magnetic field on one border is therefore an interesting 
problem because it generates a localised interface whose extension inside the bulk 
can be studied in the framework of diffusion processes. In particular the interface 
generated by two opposite magnetic fields of infinite amplitude in the critical Ising 
model is a typical example of Schramm-Locwner evolutions (SLE) ]T^ . It is generally 
restricted to the conditions of infinite strip models, and the typical size of the interface 
is scaled with the strip width. In this paper we want to focus on another particular 
geometry that can be treated exactly, where we apply a non homogeneous surface 
field on a rectangular Ising lattice. We will see that the interfacial free energy can be 
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Figure 1. Inhomogeneous field configuration on the lattice. The field is positive 
for the sites {m,n) = (1, 1) . . . (1, Ly/2) and negative on the rest of the line. 
Periodic conditions are imposed along the y-direction, and open conditions on the 
transverse direction. 



evaluated in the discrete case and that, in the thermodynamic hmit, two dominant 
and competing terms determine whether the interface stays locahsed on the surface or 
is extended inside the bulk, depending on the lattice aspect ratio. This model, even if 
the perturbation stays on the surface, shows some interesting complexity in the bulk. 
The techniques used here are based on Grassmann variables ^1 El El El El E] 
and can interestingly be extended to other configurations as well. Experimentally it 
might be important to take account of the geometry of the system like the ratio of 
typical length scales to determine the values of the critical field or temperature, as we 
are going to see in the following. 

This article is organised as follow: We summarise in section O a previous exact result 
in the discrete case of the free energy contribution due to a inhomogeneous magnetic 
field. The analytical expression is however not easy to deal with and is rather obscure 
from a physical point of view, even if it has the advantage to prevent the use of any 
cut-off, inherent to continuous models. Simple arguments for the formation of the bulk 
interface as function of field amplitude are given at zero temperature. An asymptotic 
method to obtain the thermodynamic limit of the discrete solution is presented in 
section (31 and show the simple effect of the system size ratio on the transition line. 
Monte Carlo simulations are also used to confirm the phase diagram. We then conclude 
and propose briefly the possibility of extensions of our method to other surface effects. 



2. The model 



We consider in this paper a finite 2D Ising system with a non homogeneous magnetic 
field hn located on one boundary of the system. This system is periodic along the 
y-direction, see figure^ and with open boundaries for the transverse x-direction, one 
of the latter boundaries being under a magnetic field. The Hamiltonian is simply given 

by 

H^—J {(JmnO'm+ln + CTmnCrmn+l) — ^'^^ hnCrin, (1) 

7n,n—l n—1 

with (Tmi = (JmLy+1 and con = cTL^+iri = 0. The notations are the same as in a 
previous publication |2()| where an exact expression for the free energy in the case 
of finite size and discrete lattice case was obtained using Grassmann techniques and 
Plechko method based on Grassmann operator ordering |15l I16| . These operators 
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Figure 2. Spins configurations at T = under various conditions; (a) L^/Ly > 
Cs, H < J(l + i/Ly), (b) L^/Ly > C,^, H > J(l + i/Ly) aud (c), L^/Ly < Cs, 
H > H^=4JL^/Ly. 



replace basically spin operators, and in the 2D case, they lead to a Grassmannian 
quadratic action which is exactly solvable in the Fourier space. We have found in 
particular an exact ordering of border operators associated with a general boundary 
magnetic field and then obtained, after integration over bulk degrees of freedom, a 
quadratic ID action with effective Grassmannian magnetic fields (see equation (61) 
in reference HO]), for equal system sizes = Ly. This ID action represents the free 
energy contribution from the boundary fields. The case with different sizes ^ Ly is 
easily implemented as we will see below. As an application, we have considered in j2(J) 
the interface initiated by an inhomogeneous magnetic field of configuration shown in 
figure n with hn — H ioi n ~ 1 . . . Ly/2 and /i„ — for n ~ Ly/2 + 1 . . . Ly. 
Computing the interface energy (Tint in the inhomogeneous case is equivalent to solve 
a set of Grassmannian two points correlation functions (see equation (74) in |2()j ) using 
the previous ID boundary action for the homogeneous case, which is done exactly. For 
more general configurations, with different sequences of fields {Hk'Jk), i-e. hn = Hk 
for n = Ik + \ ■ ■ ■ Ik+i, the problem is treated identically [50]. In this paper, we will 
consider the solution for the interface energy (see equation (84) in [201 and below), 
with Lx ^ Ly, and study the thermodynamic limit with fixed aspect ratio C = Lx/Ly. 
Numerically, the difficulty to compute this free energy arises from the fact that it is 
expressed as the logarithm of some argument which is an exponential small number 
in the system size, especially at low temperature. Moreover, to study in detail the 
phase transition corresponding to the propagation of the interface from where the 
magnetic field changes its sign, we need to take directly the thermodynamic limit, 
and the purpose of this paper is precisely to study how the discrete expression of the 
interface energy behaves when we take the limits Lx, Ly —>■ oo. 

The presence of an interface phase transition at zero temperature can be analysed 
with simple energetic arguments. For small values of the field H, all spins are pointing 
in the same direction, say up, because boundary negative fields are not strong enough 
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Figure 3. Boundary spin-spin correlation function (criocrir) as a function of r jLy 
and aspect ratio C, = Lx/Ly, at T = 2J, H = O.IJ and for Ly = 100. Note the 
crossover from a ID behaviour to a 2D behaviour at f ~ 1/4. 



to compete with transverse Ising couplings and reverse the corresponding spins (see 
figure The energy in this case is Eq. When increasing the field, the spins 

cl„/2+i,i • • • CLy.i will eventually reverse their sign, and the corresponding energy 
is El ^ Eq - HLy + 2J{Ly/2 + 2), which is lower than Eq if H > J(l + A/Ly) (see 
figure ISh)- Another possible configuration is when all spins <JLy/2+i,j ■ ■ ■ <^Ly.j , for 
j = 1 . . . Lx, reverse their sign (see figure In this case, the total magnetisation is 
zero, and the corresponding energy is E2 = Eq — HLy + AJL^^ which is lower than 
Eq for H > 4(J. Comparing the energies Ei and E2, we conclude that the interface 
stays on the boundary if C > 1/4 + l/Ly = Q {Ei < E2, and H > ,/(l + 4/Ly)), 
and propagates inside the bulk when C is smaller than the critical ratio value Cs and 
H larger than = 4:(J. In the latter case, where the total bulk magnetisation 
spontaneously goes from unity to zero, the transition is first order. 
In the thermodynamic limit, (^g tends to 1/4. This particular value is actually deeply 
related to the boundary condition in the y-direction. For free boundary conditions 
we would have found (^3 = 1/2 instead. The physical interpretation of this threshold 
is simplified by the study of boundary spin-spin correlation function for various C. A 
direct extension of results in [201 leads to the figure 01 For <C Ly we observe an 
exponential decay of the correlation functions, typical of a ID behaviour. If the aspect 
ratio increases, we observe a obvious crossover towards a 2D behaviour at large ^. The 
crossover is obtained in this case for C ~ 1/4. 

In particular, we would like to know in the following how the transition line Hs{T) 
for C < 1/4 behaves as the temperature is increased to near the second order phase 
transition at T^/J = 1/ tanh"^(A/2 - 1) ~ 2.27. 
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3. Analysis of the thermodynamic limit 

The Hamiltonian leads to the decomposition of the total free energy T as follows: 

T{T, H) = ToiT, H) + ^hom(T, H) + tTi„t(r, H), 

where J-q is the free energy in zero field, ^hom the additional free energy corresponding 
to an homogeneous boundary magnetic field H . The contribution tJint is the correc- 
tive term due to the change of sign of boundary magnetic field. The exact expressions 
for those terms in the discrete case can be found in |2Uj . As (Tint is the only term 
corresponding to inhomogeneous surface conditions, it contains all the physical infor- 
mations about the transition described in the previous section. The section H3.1|l starts 
with the calculation of dint in the thermodynamic limit and the evaluation of finite 
size corrections as well. This allows us to characterize the details of the transition: 
In section (|3.2|) , we obtain the expression of the transition line and the corresponding 
phase diagram, and in section H3.3|) we analyse the behaviour of this line at low tem- 
perature and close to the bulk critical point. Finally, in section (|3.4() . we summarize 
the different results and the physical interpretation. 



3.1. Expression of the interfacial free energy 



We start with the discrete expression of the boundary free energy crint taken from 
reference [50] (see equation (84) in that paper). It was obtained by computing two 
point correlation functions with a Grassmannian quadratic action expressed in Fourier 
modes. The corresponding result is the following: 



f3cr int = log 



o Ly/2-1 

— y 



(-l)'cot(0„+i/2)F(cosVi) 



(2) 



where 9g = {2iT/Ly){q + 1/2), and F is a function of < = tanh/3J and u = tanh/3i/ 
defined as 

4tu^G{x) 



F{x) 



i[l - (1 + t2)(i2 + 2tx - l)G(a;)]2 -f- 2tu^{l + x)G{x) + At^{l - a;2)G(x)2 



The function G is defined by 



G{x) 



^ Lx-l 



p=0 



{1 + fiy - 2t{l - t^){cosep + x)'^ 



(3) 



(4) 



where 9p = 2ttp/Lx. We propose in this section to simplify the expression ||2J) by 
studying the thermodynamic limit, in order to obtain the dominant terms contributing 
to the free energy. The sum inside the logarithm function jSJ has the particularity 
to behave like a Dirac distribution in the limit Ly — > 
following sum 



Indeed, we define first the 



Ly/2-l 



SiP]-T- E (-l)'cot(Vi/2)i^(cos(Vi)): 



(5) 



q=0 



for any function F. In [Appendix A| we demonstrate that in the thermodynamic limit, 
this sum is simply F{1), plus corrections which are exponentially small in the system 
size Ly. These corrections are important for finding the asymptotic behaviour of the 
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interface free energy. We are expecting a-i^t to be linear in L^. or Ly^ so the argument 
in the logarithm (j2J) should be exponentially small in or Ly. In particular, the 
Ly dependence is contained in the corrections of the Dirac distribution, and the 
dependence is contained in the term F{1) through the finite sum G, equation (0J. 
From [Appendix X] and equation (|A.7p . we can expand S[F] in the limit of large Ly 
and we obtain for the dominant part 

S[F] ^ ^^(1) - Alj2, (6) 

where ^_l„/2 is the Ly/2-th Fourier coefficient of the function F, and it will be shown 
that it goes to zero exponentially in Ly. The term F{1) is equal to 

^ ' [l-(l + t2)(t2 + 2t-l)G(l)]2/4 + 4tu2G(l)' 
and G{1) given by can be expanded as a Fourier series the following way 

Lx 

G(l) = — ^^i?,cos(fc0p), (7) 

^ p=0 fc>o 

where the Fourier coefficients Bk are defined by 

1 1 



Using 



271 Jo '^^(I+t2)2-2t(l-i2)(l+coS0)' 



TT Jo (l+t2)2 _2t(l-t2)(l+cOs6») 



1 /■^'^ d0 1 



(g) 

2tt Jo a-hcosB y/o? - ' ^ ' 

we arrive at the expression Bo = I/[(l + t2)|t2 + 2t — It is then easy to show that 
G(I) = E ^fc;^ E ) = ^0 + + ^2L. + . . . (10) 

fe>0 ^ p=0 

We expand then F{1) for La; large, when Bl^ is small, and it is sufficient to keep the 
first two terms: 

p(^. ^ 4tu^{Bo + BLj 

ill - (1 + i2)(t2 + 2t - l)(i?o + i?Lj]2 + 4t^i2(B„ + S^J ' 

^^_0^fnt^ + 2t-llBl 

16tw2 Bo ' ^ 

We have check that the coefficient B2l^, which would be of the order of , does 
not appear at this order. Bj^^ can be computed analytically in the complex plane. 
If we define z = cxp{i9), we have Bl^ = Re ^ dz/(2z7r)4z^^/(3(z), where Q{z) is a 
polynomial function 

Q{z) ^ - 2t(l - t^)z^ + 2z[{l + i2)2 _ 2t(l - ^2)] _ 2t(I - t^), (12) 

= -2t{l-t^){z-z^)(z-z+). 

The zeros z± of this function are distributed on the positive real axis with z_ = 
(1 - t)/[t{l + t)], z+ = l/z_. z_ is less than 1 (or z_ < z+) for t > ic = \/2 - 1, 
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in the low temperature regime (therefore the quantity + 2t — 1 is always positive) . 
The value of B^^ in this region is then equal to 

" (l+i2)(t2 + 2t-l) (i(l+t)) ■ ^^^^ 

We obtain therefore, for the dominant part of the distribution S: 

These corrective terms are all negative, since 1 — S[F] should be positive so that the 
logarithm in (j^J is always defined, and the interfacial free energy can now be written 
as 

-Paint = log(l-5[F]), 

, f {l + t^){t^ + 2t-l) f l-t ^ \ 

^ y- — hi^^ — - wtt) ) + ^-"/^ ) ■ (1^) 

The Fourier coefficient is evaluated for large Ly with the function 0. We 

expect this coefhcient to be exponentially small with Ly, with some corrections to this 
coefficient which are also small in Lx- Therefore the dominant term can be obtained 
by taking the limit L^; = oo in Q). Using the formula © we obtain in this limit 

G{cose)^ J (16) 

R{cos9) = [(f + i2)2 + 2t{l - t2)(i - cos6')][(f + t^f - 2t{l - + cos 6')]. 
The function F can then be rewritten, after some algebra, as 

8f7/2 

F{cos9)^ — , (17) 

Atu^{l + cose) + {l + t^){l~2t cos e-t^) + ^/R{cose) 

Now ^Lj,/2 can be expressed by mean of a complex integration along the unit circle 
AlJ2 = Rei ^ d0^^(cos 0) exp (^z^^^^ 

The value of this integral depends on the poles of the function F. Setting X = 
{z^ + l)/2z, we can write 

Fix) = = !^ L fig-) 

^ ' P{X) + y^RiX) P{XY - R{X) ■ ^ 

The polynomial P{XY — R{X) is a second order polynomial in X, and has two zeros 
which are 

1 2t^-2tu^-u\l~t^) 
^'-2 t[l~u^W-u^) ' (20) 



Xi is not a pole of the function F since it is not a zero of P{X) + ^ R{X). Indeed, 
we have P(-f ) + ^/R{-1) = {\+t^){\ + 2t-t^ + \l + 2t-t^\) > 0, because I + 2i - 
is always positive in the interval < t < 1. Then only Xq is pole of F, and in the 
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<M>=0 




0,5 1 1,5 2 2,5 3 3,5 4 



Figure 4. Phase diagram for the system at (J = 0.2. The plain Une represents the 
first order transition given by equation 1231 . The dashed hne is the bulk 2"^* order 
phase transition. The snapshots are extracted from Monte Carlo simulations in 
the different regimes, for = 40 and Ly = 200. 



complex plane this gives two solutions Z± of the equation — 2Xqz +1 = 0. If Z_ is 
the solution whieh is inside the unit circle {\Xq\ > 1), the value of is given by 



p{Xo) - VRi^) 



{p{xY~Q{x)y{x,)^l 



1 



(21) 



In some cases the zeros can be on the unit circle {\Xq\ < 1), but we assume this latter 
situation is not physical since this would mean that the thermodynamic limit is never 
obtained. 

Finally the free energy for the interface in the large system size limit can be written 
as the logarithm of two dominant terms 

+ C,(t,«)|Z_|^='/2 ^ (22) 



/3(Tint ^ log 



+ + 2t-i) 



1 - 1 

i(l + t) 



3.2. Line transition and phase diagram 

Two different regimes can be identified from equation (|22|l . depending on whether 
the first term in the logarithm is greater or smaller than the second one. In the first 
case, the free energy is proportional to Lx and a coefficient which does not depend on 
magnetic field, only on t. At zero temperature, it is easy to show that CTint — 4:J(^Ly, 
which corresponds to the energy E2 computed above in section This means that 
this term represents a situation where it is energetically favourable for the interface 
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to spread inside the bulk. On the other hand, the second term in the logarithm 
of equation (|22|l represents a configuration where the interface is localised on the 
boundary. Between these two regimes, we can define a line of transition which is a 
priori first order: The first term does not depend on the magnetic field, therefore the 
free energy has in general a cusp as a function of u) by making the two exponential 
terms in (j22|l equal in magnitude. 

We then obtain the transition line equation in the {t, M)-plane which is simply a quartic 
polynomial in u: 



2t I 



w(4C)] m"^ + (1 + t^) (l - 2to(4C) 



0, 



(23) 



v{A() = cosh 



4C In 



1-t 



On figure 21 is shown the phase diagram for the system at C = 0.2, with snapshots of 
typical configurations at low and large fields/temperatures, obtained by Monte Carlo 
numerical simulations. 



3.3. Behaviour around T = and T ^ 



Near zero temperature, we can expand the previous relation with f = (1 — e)/(l + e), 
e = exp(— 2/r) <C 1. We obtain at lowest order in e and for 4^ < 1: 

(1 + 2e^^)u^ -2u^ + l- 2e^'^ = 0, (24) 

which gives, = 1 or = (1 — 2e"'^)/(l + 2e'^''). The non trivial solution gives the 
point Hg — 4C,J as expected at zero temperature from preliminary study. Moreover, 
if the transition line ends at the point w = 0, this is equivalent to Xq = u = 1, or 
1 — t = t{l + t), which gives t = tc. The line ends therefore at the second order 
transition point. In this case Z+ = Z_ = 1, which is the transition value between 
an exponential behaviour and oscillating one in the logarithm arguments (j^J. This 
basically suggests that the interface free energy is no more an extensive function of the 
system size, and does not contribute to the thermodynamic behaviour of the system. 
In the case when 4^ > 1, the discriminant A of the equation (|23() can be written as 



A = At^{l - f 



4t2(l - f 



{i + ey + 2t{i-e) 

2^(1-^2) 

{i + ef + 2t{i-e) 



(l+t2)2-2t(l-P) 



2t(l-t2) 



{v-v{l)). 



2t(l-t2) 

Expanding v{AC) near the threshold value 4^s = 1, we obtain 



(25) 



w(4C)~t;(l) + (4C-l)ln 



l-t 



sinh In 



_<(l + t) 

and the discriminant can be expanded as 



1 - t 



t{l + t) 



A ~ 8t^(l - t^yin 



1 - t 



t{l + t) 



sinh In 



1 - t 



(1-4C). 



(26) 



t{l + t)_ 

The discriminant A is negative when > 1/4 and equation ()23(l has no real solution, 
therefore the wetting transition no longer exists in this regime. 

Near zero temperature and for <^ < 1/4, we can expand more precisely equation H23|l . 
for small parameter e, this will give locally the behaviour of Hs{T) as function of 
temperature. We obtain for the quantity v: 

v = ie-^^(l- 8Ce + e^'^ + ...). 
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This has to be equal to the expansion of Xq, which is given, at lowest order, by 

Xo = le-"il + e^" + ...). 

By comparing the two previous expressions, we obtain Hs{T) ~ 4^ — 4^Texp(— 2/T). 
Figure^ shows that, as expected, the curve Hs{T) is flat near zero temperature due 
to exponentially small thermal excitations. Near the critical temperature tc, a simple 
analysis gives Hs{T) cx y't — tc- 



3.4- Summary of the results 

To summarize the previous calculations, we have found that the free energy can 
be expressed as the logarithm of two contributions which arc exponentially small in, 
respectively, the system sizes and Ly, see equation (|22|l . This result is obtained after 
performing an asymptotic analysis of the Fourier sum S, equation |(5j, that appears 
in the logarithm of equation Q. This sum behaves like a Dirac distribution in the 
thermodynamic limit (see |Appcndix A| ) and its value in this limit makes the logarithm 
in |5J) singular. To remove the singularity we analysed the finite size dependent 
corrective terms since we expect (Tint to diverge linearly with the system size. The two 
main contributions inside the logarithm essentially come from two relevant Fourier 
amplitudes that tend to zero exponentially with and Ly. The other contributions 
are much smaller and do not contribute to the free energy. One term corresponds 
to the interface localized on the boundary and the other to the interface extended 
across the bulk. The relative amplitude between these two terms is controlled by 
the aspect ratio C, and the transition line between the two regimes is expressed by a 
simple quartic equation l|23() . A study of its discriminant H26|) leads to the existence of 
a first order transition line when C < 1/4. For ( > 1/4 no real solution exists, and the 
interface is always localized on the boundary since only one of the two contributions 
inside the logarithm is dominant at all temperatures and magnetic fields. 



4. Conclusion 



In this article we obtained an exact description of a first order phase transition 
induced by a simple inhomogeneous boundary magnetic field. The use of Grassmann 
techniques allows for an exact calculation of the interfacial free energy in the discrete 
case, which is suitable to study then the thermodynamic limit by asymptotic methods 
presented in this paper. This approach allows us to control the way the thermodynamic 
limit is taken, and has the advantage that no cut-off parameter is required for the 
continuum limit. In particular, we have demonstrate that it is straightforward to take 
the thermodynamic limit exactly for a given geometry. This leads to a surprisingly 
simple equation of the transition line H28(l in the {H, T) or (u, t) planes, and the 
corresponding critical behaviour close to the bulk critical point (0, Tc). In the context 
of wetting transitions this is a non trivial extension of previous results. This line 
disappears for ^ > 1/4 as the solutions move to the complex plane. A numerical study 
of bulk correlation functions at the precise value C = 1/4 might show the dynamical 
instability of the interface on the vanishing transition line. The infinite strip models 
for inhomogeneous surface field |Z| El E| might not capture this feature since the 
interface is not sensitive to the system aspect ratio. 

Grassmann techniques, in complement to conformal theory |21l I22j and transfer 
matrix methods, can be considered as an interesting optional way to solve boundary 
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problems or wetting transitions. Extensions of the method presented here and in j2(J) 
might be useful to study other kind of wetting transitions. In particular it might 
be applied to models of defects other than a surface field such as a line of weaker 
or stronger couplings |23[ , which has been studied in the framework of transfer 
matrix methods in detail, and where striking similarities with other physical domains 
like electrostatics have been suggested. 

Appendix A. 

We would like to compute the following sum in the large L (even) limit (Ly is replaced 
by L here for generality): 

L/2-1 



^[^] = 7 E (-l)'cot(0,+ i/2)F cosVi 



9=0 



(A.l) 



with ^q+i = (27r)/L(g + ^), and F any function of the variable cos^^^i . We know 
that for any constant F, S[F] = F (see [Appendix B| . For commodity, we can extend 
the sum from q = 0..L/2 — 1 to q = 0..L — 1 hy writing 

L/2-1 



S\F] 



1 2 
2L 



E (-l)'cot 



q=Q 



F 



1 9 



cot 



q=L/2 



F 



, 2n , 

cos[ — [L-q-l + - 



It is then straightforward to see that, after rearranging the different terms 

L-l 



7=0 



T\<1 + 1 



F 



cos [ - [ q + i 



(A.2) 



Next we express the function F as a Fourier series: 



F(cos 6*) = E cos(p6i). (A.3) 

p=0 

Computing S[F] is equivalent to obtain explicitly the value of every term S[cos{p9)]. 
For p = 1, we have 



cot I 



cos 9„ 



1 + cos I 



sin 9„,i 

9 ' 2 



sinfl 



COSt/^+l 

— sin d„. 



COS 0„ I 1 + cos^ 9„^ 
sin 9„,i 

2 



9+5 



which implies that 

5[cos(.)] = 5[l]-iE(-irsin0,+ i. 



L-l 



9=0 



For the second term, p = 2, we can show that 

1 + cos6'„ , 1 
cos(20„+i) = . ^ 

^ 9T^ 2 ' cin fl 



1 

cot I - 



cos(26'„, i). 



cos6'^_i_i + 1 



sm 



cos 6l^_^i - sin(26l^+i) - sin 61^+1 , 



'9+ 



(A.4) 
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and for the third term 



/I \ l + cos6'„,i 

cot -9^,1 cos(36i„, i) = — — — ^cos(36'„, i), 

cos 6„.i + 1 

By recursion, it is easy to show that 

5[cos(p.)] = 5[cos([p - 1].)] -Tp- Tp_i, (A.5) 

FinaUy, after rearranging the previous relations, we obtain the foUowing expression 

5[cos(p.)] =S[l]-Tp-2(rp_i+Tp_2 + ... + Ti), 5[1] = 1, (A.6) 

5'[cos(.)] =5[1]-Ti. 

The quantities Tp can be easily computed, and we obtain Tp = except for 
p = i(i + k), k > 0, where Tp = (—1)''. After some combinatorial, the distribution 
S[F] is equal to 

oo oc L — 1 

S[F]^F{l)-Y,AL^2+kL-2Y^ AL/,+2kL+k'. (A.7) 

k=0 fc=Ofc' = l 

In the limit where L is infinite, S[F] is the Dirac distribution S[F] = F{1) since all 
the Fourier coefficients tend to zero. 



Appendix B. 



In this section, we want to prove the following equality 

L-l 



q=0 



TT , 1 , 



= 1, 



(B.l) 



for any value of L. Let assume that L is even, the proof for L odd is equivalent. We 
can notice that 

L-l 



d_ 



log 



q=Q 



sm 



L 



Separating in the sum the odd and even integers q, we obtain easily 

L/2-1 



S[l] 



ld_ 
T'd'z 



log 



n 



sin[7r(2g+ l/2)/L- 



L sin[7r(2g + 3/2)/L + z] 



(B.2) 



(B.3) 



The 2 products are evaluated after expressing the sine function as exponential terms, 
and using by identification the equahty X^l"^ — 1 = Ilgio^^ ~ cxp(4i7r(7/L)]. We 
then obtain the simple result 

sin(zL/2 + 7r/4) 



1 d 

5[l] = --log 
L oz 



sin(zL/2 + 37r/4) 



1 



[cot(7r/4) - cot(37r/4)] 



z=0 
1. 



(B.4) 
(B.5) 



By extension, S[C] — C for any constant C. 
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